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Abstract. We consider the integrable multicomponent coherently coupled nonlinear 
Schrodinger (CCNLS) equations describing simultaneous propagation of multiple 
^vq . fields in Kerr type nonlinear media. The correct bilinear equations of m-CCNLS 

OA I equations are obtained by using a non-standard type of Hirota's bilinearization 

method and the more general bright one solitons with single hump and double hump 
^^ , profiles including special flat-top profiles are obtained. The solitons are classified 

Ph ' as coherently coupled solitons and incoherently coupled solitons depending upon the 

r^ . presence and absence of coherent nonlinearity arising due to the existence of the co- 

propagating modes/components. Further, the more general two-soliton solutions are 
obtained by using this non-standard bilinearization approach and various fascinating 
collision dynamics are pointed out. Particularly, we demonstrate that the collision 
among coherently coupled soliton and incoherently coupled soliton displays a non- 
-Y-\ ' trivial collision behaviour in which the former always undergoes energy switching 

f^ , accompanied by an amplitude dependent phase-shift and change in the relative 

^O ' separation distance, leaving the latter unaltered. But the collision between coherently 

^f. . coupled solitons alone is found to be standard elastic collision. Our study also reveals 

C^ I the important fact that the collision between incoherently coupled solitons arising 

^"^ ■ in the m-CCNLS system with ttt, = 2 is always elastic, whereas for tti > 2 the 

collision becomes intricate and for this case the m-CCNLS system exhibits interesting 

energy sharing collision of solitons characterized by intensity redistribution, amplitude 

• "^ ■ dependent phase-shift and change in relative separation distance which is similar to 

r> I that of the multicomponent Manakov soliton collisions. This suggests that the m- 

C^ • CCNLS system can also be a suitable candidate for soliton collision based optical 

computing in addition to the Manakov system. 
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1. Introduction 

Multicomponent solitons/solitary waves have attracted considerable attention in the 
field of nonlinear science as they display a rich variety of propagation and collision 
properties which are not possible in their single component counterparts [IHT8]. Such 
solitons appear in different areas of science like nonlinear optics [IHS], Bose-Einstein 
condensates [121, bio-physics [20], plasma physics |3], etc. Here, our main focus is on 
nonlinear optics. In this context, multicomponent temporal solitons can be formed 
when an optical pulse propagating through a multimode fiber due to a delicate balance 
between dispersion and Kerr nonlinearity [Ij. Multicomponent spatial solitons are self- 
trapped optical beams that result from an interplay between diffraction and nonlinearity 

Mathematically, the propagation and collision properties of multicomponent 
solitons/solitary waves arising in the field of nonlinear optics can be well described 
within the framework of multicomponent nonlinear Schrodinger (NLS) type equations 
[I112]. Especially, the short pulse propagation in polarization maintaining multimode 
birefringent fiber is governed by a set of multicomponent incoherently coupled NLS 
(ICNLS) equations [21]. Similar set of ICNLS equations also arises in the context of 
partially incoherent beam propagation in Kerr type nonlinear media [TT]. These ICNLS 
equations involve the nonlinear couplings due to self-phase modulation (SPM) and cross- 
phase modulation (XPM) and depend only on the local intensities of the co-propagating 
fields, but insensitive to their phases [1]. 

In general cases, like pico-second pulse propagation in non-ideal low birefringent 
multimode fibers or beam propagation in weakly anisotropic Kerr type nonlinear media, 
the coherent effects due to the interaction of co-propagating fields should also be 
considered [HET]. To be specific, the propagation of coherently coupled orthogonally 
polarized waveguide modes in Kerr type nonlinear medium is governed by the following 
2-component coherently coupled NLS (CCNLS) type equations [T]: 

%,2 + Sqi,tt - mi + (kiP + cr\q2\'^)qi + Xq^ql = 0, 

iq2,z + 5q2,tt + /ig2 + (cr|giP + \q2\'^)q2 + MIqI = o, (i) 

where gi and g2 are slowly varying complex amplitudes in each polarization mode, z and 
t are the propagation direction and transverse direction, respectively, fi is the degree 
of birefringence. Here, the incoherent and coherent couplings are represented by the 
parameters a and A, respectively. The above equation ([1]) also arises in the context of 
beam propagation in isotropic Kerr type nonlinear gyrotropic medium [22] • In equation 
([1]), the terms g^Q'i and qfq2 correspond to four wave mixing (FWM) process which arise 
due to the coherent coupling between the co-propagating fields. 

Generally, these CCNLS equations and also the ICNLS equations are non- 
integrable. However, these become integrable for specific choices which are of physical 
significance [T ^[25H25] . In recent years, much attention has been paid to the integrable 
and non-integrable coupled ICNLS equations and many interesting phenomena have 
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been explored [1]. To be specific, the integrable m-component Manakov type equations, 
with arbitrary m, are well studied and it has been pointed out that these equations 
support bright optical solitons which undergo fascinating energy sharing collisions that 
have immediate technological applications in the context of collision based optical 
computing [ll[TTl[T2l|271|28] and also in soliton amplification ^6]. Very recently, CCNLS 
equations have also been started to receive renewed attention due to their rich structure 
[SIII31I251ES] ■ Particularly, a set of physically interesting integrable 2-component CCNLS 
equations related to ([1]) is 

iqiz + qm + i{\qi? + '^\q2?)qi - iqhl = o, 

iq2z + q2tt + 7(2|giP + \q2?)q2 - iqlql = o, (2) 

where 7 is the coupling coefficient. In gyrotropic nonlinear medium, the above equation 
(2) can be obtained for the following choice of the susceptibility tensor x^^\ with its 
components satisfying the relation XxLx = Xxlyy = Xxyxy = Xxyyx [22] • Equation (2) also 
describes the propagation of two optical pulses in an isotropic nonlinear Kerr medium 
when the components Xxxyy-, Xxyxy and Xxyyx of the susceptibility tensor x'-^'' can be 

expreSSeO. as Xxxyy ~r Xxyxy ~— ^Xxyyx j-l-OJ . 

Apart from the 2-component CCNLS equations, ?7i-component CCNLS equations 
with m > 2 are also of special physical interest and have been derived under different 
physical contexts. Particularly, the spatial evolution of mutually guided four wave 
mixing states in x^^'' medium is governed by 3-component CCNLS type equation [3D] . It 
has also been shown that the co-propagation of two optical pulses in birefringent fiber 
can be described by 4-component CCNLS equations [31j. In ref. [Tl|, the dynamics of 
spinor Bose-Einstein condensates has been investigated by considering a set of integrable 
3-coupled CCNLS type equations [26] and novel polar and ferromagnetic solitons have 
been reported. Now it is of interest to investigate the integrable multicomponent CCNLS 
equations which are closely associated with the near-integrable or non-integrable systems 
appearing in nonlinear optics. 

Being motivated by these reasons, we consider the following integrable m- 
component generalization of ([2]) describing the simultaneous propagation of m-optical 
fields in Kerr type nonlinear media. 

(m \ m 

l^il^ + 2 Yl l^«n?i-7 Yl ^i^^i=0, j = 1,2,3, ...,m. (3) 
i=Ujtj J 1=1,1^3 

The above system has been studied in ref. [25] by applying the Hirota's direct method 
but trivial soliton solutions with less number of parameters only have been reported due 
to the restricted bilinearization of ([3]). It also should be noticed that the information 
regarding the coherent and incoherent contributions from the co-propagating fields are 
lost completely if the two-soliton solution is constructed by a linear superposition as 
pointed out in refs. [T31I25] . So it is of importance to obtain correct bilinear equations of 
system ([3]) which will result in more general soliton solutions with interesting properties. 
In ref. [9], Kanna et al have considered a 2-component integrable model which can be 
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reduced from ([2]) by redefining gi as iqi and reported novel solitons with variable profiles 
and classify them as degenerate (solitons with same intensity in both components) and 
non-degenerate (solitons with different intensity in two components) solitons. But study 
on the present system (3) suggests that a broader classification of solitons of the general 
m-component system ([3]), with arbitrary m, can be made based on the presence and 
absence of the coherent nonlinearity, and the degenerate and non-degenerate solitons 
reported in [H] appear as their sub-cases, which will be discussed in the following sections. 

The aim of the present work is three-folded. First, to obtain the correct bilinear 
equations of system ^ and to construct exact one- and two-bright soliton solutions 
of ([3]). Next, to analyse the collision dynamics of solitons in the 2-component and 
3-component CCNLS equations and to bring out their salient features. Finally, to 
generalize the results to the ?7i-component case. 

This paper is set out as follows. The correct bilinear equations of ([3]) are obtained 
and the solitons are classified in a systematic way in section 2. The bright one- and two- 
soliton solutions of 2-component and 3-component CCNLS equations are obtained in 
section 3 and in section 4, respectively. Then the results are generalized to arbitrary m- 
component case in section 5. Section 6 deals with the collision dynamics of the solitons. 
Final section is allotted for conclusion. 

2. Non-standard bilinearization and classification of solitons of integrable 
multicomponent CCNLS system ([3]) 

Hirota's direct bilinearization method is one of the powerful techniques to construct 
soliton solutions of integrable nonlinear evolution equations [22]. In this section, we 
construct the bilinear equations of the CCNLS system ([3]) by applying the Hirota's direct 
method [32] . A new type of bilinearization procedure has been developed by introducing 
an auxiliary function for the Sasa-Satsuma higher order nonlinear Schrodinger equations 
in ref. [33] by Gilson et al. By adopting this technique, here we obtain correct bilinear 
equations of system ([3]) resulting in more general bright soliton solutions which display 
the effects of both intensity and phase dependent nonlinearities. By performing the 
bilinearizing transformation 

g{3) 

qj = —, j = l,2, ...,m, (4) 

to equation ([3]) with the introduction of an auxiliary function s, we arrive at the following 
set of bilinear equations: 

j = l,2,...,m, (5a) 

{5b) 

(5c) 
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where Di = iDz + Df and D2 = Df. Here g'^^'> and / are complex and real functions, 
respectively, * denotes the complex conjugate, D^ and Dt are the well known Hirota's 
D-operators [32] which are defined as 

The above set of equations (5) can be solved by introducing the following power series 
expansions for g^^\ f, and s 

9 = X9i^ + X^93^ + ■■■, j = l,2,...,m, (6a) 

/ = 1 + XV2 + xV4 + . . . , s = x's2 + x's4 + . . . , (66) 

where x is the formal power series expansion parameter. The resulting set of linear 

partial differential equations after collecting the terms with the same powers in x, can 

be solved recursively to obtain the forms of g^^\ /, and s. 

It can be inferred from the above bilinear equations (5) that when the auxiliary 
function "s" becomes zero, the contribution from the coherent coupling vanishes and the 
above bilinear equations reduce to that of integrable ?Ti-component Manakov system. 

•m 

We notice that for s = equation fl5cj) results in the condition /^(f? )^ = 0, 

which ultimately restricts the energy sharing of a given field/soliton among all its 
components. In the following, we obtain explicit conditions on soliton parameters for 
which "s" becomes zero and we refer the soliton arising for this choice, s = 0, as 
incoherently coupled soliton (ICS) as the contribution from the coherent nonlinearity 
is absent. The ICS results due to the interplay between dispersion/diffraction and 
the nonlinearity arising due to SPM and XPM effects. However the general bilinear 
equations (5) with non- vanishing auxiliary function "s" involve the effect of coherent 
coupling also. Hence, we designate the soliton resulting for the general choice, s 7^ 0, as 
coherently coupled soliton (CCS). These CCSs are formed due to the contribution from 
the dispersion/diffraction and the combined nonlinear effect resulting from SPM, XPM 
and four wave mixing process. 

3. Bright soliton solutions of 2— component CCNLS equations 

In this section, the bright one- and two-soliton solutions of the 2-component CCNLS 
equations ([2]) are obtained by applying the non-standard type of Hirota's bilinearization 
method explained in the previous section. We present the results for the m = 2 and 
m = 3 cases explicitly in order to emphasize the additional features of the m = 3 case 
{m represents the components). Here onwards we designate the m-component A^-soliton 
solution as {m, N) soliton solution for convenience. 

3.1. Bright (2,1) soliton solution 

To obtain bright one-soliton solution of system (2), we restrict the power series expansion 
(6) as c/(^) = X9i^ + X^9?^ J = 1, 2, / = 1 + XV2 + xV4, s = x^S2. By substituting 



where 
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these series expansions into (5) and after recursively solving the equations resulting at 
like powers of Xi "we obtain the following one-soliton solution. 

^J' ~ 1 + Qm+ril+Ri + e2r;i+2r,*+eii ' ^ " 1' 2, (7a) 

4(fci + A;*)4 ' ""'' {h + kl) ' ^'"^ 

The auxiliary function s is found to be 

s = {{a^^f + {af^f)e'^K {Id) 

Here, r/i = fci(t + iA;i2;), fci = kui + iku, and a]^^ s are complex parameters. Throughout 
this paper, the real and imaginary parts of a parameter are represented by the subscripts 
R and I, respectively, 
(i) Bright (2,1) ICS: 

The (2,1) ICS results for the vanishing auxiliary function [s = 0). We find from 
fl7^ that the condition for s to be zero is {a[ )^ + {a[ ')'^ = 0. This (2,1) ICS always 
exhibits the standard "sech" type profile and can be expressed as 

q, = A, sech {ri^R + R^/2) e'"" , j = 1, 2, (8) 

where Aj = ^^^=, Ri = In l^),ViR = kiR{t-2kuz), and r/u = kiit+{klj^-klj)z. 

V27l"i P ^ ^ 

For this case, either a[ = ±ia[ or a]^ = ±iai and correspondingly the solitons in 
gi and ^2 components are related as g2 = ±^Q'i or gi = ±iq2- Ultimately, the intensity 
profiles of these ICSs are same in both the components (that is, |gip = |g2p)- One can 
also refer these equal intensity solitons in both components as degenerate (2,1) ICSs 
and are characterized by two complex parameters ki and a\ (or a[ ). These solitons 
behave as standard NLS solitons during propagation. The amplitude of soliton in the 
g^-th component is Aj. The velocity and central position of soliton in both components 
are 2kii and -^^^ respectively, and can be tuned by altering either ki or a\ {or a[ ). 
Such an incoherently coupled soliton is depicted in figure [1] for the parameters 7 = 2, 
ki = 1 — i, a\ = 1 + i, and a\ = I — i. 
(ii) Bright (2,1) CCS: 

The (2,1) CCS solution can be obtained for non-zero auxiliary function {s 7^ 0). 
We require {a[ )^ + {a\ )^ 7^ 0, for non-vanishing s and the corresponding (2, 1) CCS 
solution is found to be 

f cosjPj) cosh{Q)+ism{Pj) smh{Q) \ ,^^^ , _ 1 9 (Q^ 

V 4cosh''(Q) + L J 

where A, = e ^ , P, = '-^^, I, = ln(aS^O, J = 1,2, Q = r],R + ^, 

L = e(^i-^)-2, T]iji = kiR{t-2kiiz), and r^u = kut + {kl^-klj)z. The (2,1) CCS can 
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Figure 1. Degenerate (2,1) incoherently coupled soliton. 

exhibit both equal and non-equal intensities in both components, which may be referred 
as degenerate and non-degenerate (2,1) CCSs, respectively. Thus the degenerate and 
non-degenerate solitons obtained by one of the authors and co-workers in ref. [H] can 
be deduced as sub-cases of coherently coupled solitons of system ([2]) discussed here. 
Generally, these CCSs admit double hump profiles. The existence of coherent coupling 
is reflected by such kind of distinct profiles. One can also obtain perfect "sech" type 
soliton profile when the parameters are chosen suitably. This can be achieved for the 



condition q\ 



(1)* (2) 



a\ — a^ 



(1)^,(2)* 



a 



0, which makes Pj = L = 0, j = 1, 2. Here Aj 



represents the amplitude (peak value of the envelope) of soliton in g^-th component and 
2kij and ^^ are the velocity and the central position of the soliton in both components, 

respectively. These (2,1) CCSs are characterized by three complex parameters /ci, al 
and a\ . The non-degenerate (2,1) CCS having double hump profile in gi component 
and special flat-top profile in g2 component is depicted in figure [2] for the parameters 
7 = 2, /ci = 1—i, ai = 1, and ai = lAi. Such flat-top type solitons have been reported 
in non-integrable complex Ginzburg-Landau equations [3]. Thus from our above analysis 
we observe that one can switch from coherently coupled soliton to incoherently coupled 
soliton and vice- versa by tuning the polarization parameters (a'^-'^'s) suitably. 



3.2. Bright (2,2) soliton solution 



The bright two-soliton solution of system ([2]) is obtained by restricting the power series 
expansion (6) as g^^'^ = XQi^ + X^9? + X^95^ + X^^{^\ J = 1, 2, / = 1 + ^^2 + X^U + 
X^fe + X^/s) s = x^S2 + x'^Si + x^Sq. Then, by solving the resultant set of linear partial 
differential equations, we get the bright two-soliton solution as 



?i 



9 



U) 



J = 1,2, 



;ioa) 
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Figure 2. Non-degenerate (2,1) coherently coupled soliton with double hump and 
flat-top profiles. 



where 



.(i) 



= «(^')e''i + a^^^e''" + g^^'i+^'i+^n + (P^m+vl+^i-^ + g2»)2+'?i*+4i' _^ g2»?2-i-»)|-h4i' 



+ e 



+ e^ 



,m- 



i+»?i*+r?2+<5<'' _^ gr?2+r?2*+r?i+4'' ^ g2r?i+2r?*+»72+MH' _^ g2r;i+2r?2*-H»?2-HM^i' 
_^ g2»72-H2r?*-H»7i+4i' _^ g2»72-(-2»75+r,i+4'2' _^ g2r?i+r?*+r?2+»72+Mi'' 

g2»72+»?5+»7i+^i+4'' + g2r?i+2»7*+2»72+r?2*-h</.p' _^ g2»7i+2r?2+2r?2*+»?i*+4'' 1 = 12 



/ _ ]_ _|_ gm+^iH-^i _|_ g^i+'?2+'5o _|_ g'?2+^i+'5o _|_ g'?2+'?2+-R2 _|_ g2»7i+2r?*+eii 



+ e- 



,2,7 



(106) 



\\+2ri2+ei2 _|_ g2r;2+2»7j+e2i i g2r;2+2r;2+e22 i 2»7i-(-r;J;+»72+ri , 2r;J;+»7i+r;2+-rj* 
_j_ ^2r)2+r)l+r)l+T2 _|_ g2r?5+»7i+r?2+-r2* _|_ g»?i+»?i+»?2+r?2+i?3 _|_ g2»?i+2j7*+r?2+»?2 -1-611 
_|_ g2j7i+2r?5+»72-Fr?*-Fei2 j^ ^2ri2+2ri*+'ni+ri^+e2i ^ g2r?2+2r?5 +771 -Fry* +6*22 _|_ g2(r?i-|-r?*-Fr?2+»?5)+i?4 (IQc) 
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and the auxiliary function s is given by 

I g'')i+»yi+2»;2+Aii I prii+ril+2ri2+\i2 i g»72+''yi+2r;i+A2i i g-'?2+?72+2»?i+A22 
_|_ g2r?i+2r?J+2r?2+Ai _|_ g2-^i+2r?2+2r7*+A2 _|_ g2j7i+r?J+2r?2+-»?|+A3 (lOc/) 

Here, r/; = ki{t + ikiz), I = 1,2. Various other quantities appearing in the above equation 
(10) can be obtained from the Appendix by substituting m = 2. The above general two- 
soliton solution is characterized by six complex parameters ki, A;2, a[ , eta , a[' and 

(2) 

ay. 

4. Bright soliton solutions of 3-component CCNLS equations 

We obtain the exact bright one- and two-soliton solutions of three component CCNLS 
equations in this section by applying the non-standard Hirota's bilinearization method 
described in section 2. 

4.I. Bright (3,1) soliton solution 

The bright (3,1) soliton solution of CCNLS equation ^ with m = 3, can be obtained 
by terminating the power series expansion (6) as g'^^'' = X9i + X^9i i J = l5 2,3, 
/ = 1 + x^/2 + x'^/4; -5 = x^'S2 and by solving the set of partial differential equations 
arising at like powers of x- Then the (3,1) soliton solution can be written as 

Qi = , « , p t; — rrr^, , 7 = 1,2,3, (11a) 



where 



7a 



?*E(«!' 



e4i' = '=1 , e^^ = -^ii— . J = 1,2,3, (116) 

2{ki + kiy ' {ki + kl)' -^ ' ' ' ^ ' 



r 



t(-?r^ ^t\-?\' 



e"^" = r, r-Ti — , ^^ii = —r, r^—- (He) 

4{ki + kl)^ {k, + kl) ^ ' 

The auxiliary function s is found to be 
3 
s= ^(a;^))V^ {lid) 

i=i 
The above bright one-soliton solution can also be classified into the following ICS and 
CCS as in the previous section. 

(i) Bright (3,1) ICS: 

The (3,1) ICS appears for the choice {a\ Y + (c^i Y + (o^i Y = 0, and in its exphcit 
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qj = A, sech Ur + ^^ e^"", j = 1, 2, 3, (12) 



where Aj = a[^^^l^, Ri = In (^), ViR = hnit - 2kuz), and rm = kut + {kjj^ - 

k\i)z. Unlike in the two-component case, here the degenerate ICS is not at all possible. 
These (3,1) ICSs admit either non-degenerate (completely different intensity profiles in 
all the three components) or partially degenerate (same intensity profiles in any two 
of the components) soliton profiles. These solitons are equivalent to the 3-component 
Manakov type solitons [TU] and are characterized by three arbitrary complex parameters. 
Such non-degenerate type ICS arising for the parametric choice fci = 1 — i, 7 = 2, 

(1) _,/^ ^.(2)__/7T _, ^(3) 



a\ ' = y/2, a\ = y/3, and a\ = y/Ei is shown in figure [31 
(ii) Bright (3,1) CCS: 

The (3,1) CCSs appear for the choice {a[ )^ + {a\ )^ + {a\ )^ 7^ 0. The exact form of 
this soliton can be obtained by rewriting the one-soliton solution (I flap as 

„ _ 9 /I /^ cos(Pj) cosh(Q) + z sin(Pj) sinh(Q) ^^ ■ ^^ ,• _ 1 9 q (^'i\ 

where Aj = e 2 , Pj = '-^^, /, = Ha['>), j = 1,2,3, Q = rjiR + '-f, 

L = e(^i-^) - 2, r]iR = kiR^t - 2kuz), and riu = kut + [k]^ - klj)z. These (3,1) CCSs 
admit both single hump and double hump profiles. In fact, we obtain perfect 'sech' type 

(1) (2) (3) 

(3,1) CCSs for specific choice of parameters satisfying the relation ^^t^ = -%rj = -^tW, 

a\ al a\ 



which ultimately makes Pj = L = 0, j = 1,2,3, in the above equation flT3|) . The 
CCSs can also have degenerate intensity profiles in addition to non-degenerate profiles. 
A typical degenerate and also a non-degenerate type CCSs are shown in figure H] and 
figure Ofor the parametric choices ki = 1 — i, 'y = 2, a[ = 2, a\ = 2, and a[ = 2 
and ki = 1 — i, 'y = 2, al = 0.25, a\ = —0.71, and al = 1.25z, respectively. 

4-2. Bright (3,2) soliton solution 

The bright (3,2) soliton solution of system ([3]) with m = 3, can be obtained as in the 
two-component case by restricting the power series expansion (6) and by allowing j to 
run from 1 to 3. The two-soliton solution is found to be 
qU) 
13 = -J^ J = 1,2,3, (14a) 



W) O»„_L„*_l_^0) O^„_1_.„*_la0) O„„_1_^*_lx0) 



where 

_^ ^Vi+nl+m+A'^ _^ g»?2+r?2*+r?i+4'' _^ g2r?i+2r?*+772+A'K' _^ ^2l^l+2r)l+l^2+^l''^^ 

j^ g2,?2+2r,J+77i+4'i' _^ ^21^2+21^i+m+^l'i^ j^ g2r,i+r,J+r?2+»7|+Mi'' 

_^ g2,72 +772 +'71+^1+4'' + g2r,i+2,7*+2,72+r7*+(/.<^') _^ ^2,7i+2r72+2r7*+>7i*+4-'' 7 = 123 (146) 

f — ]_ _j_ ^m+vl+R-i _j_ ^m+'nl+^o _|_ ^n2+vl+st-, _j_ ^n2+vl+R2 _j_ g2?7i+2r7j+eii 

_|_ g277i+2r7j+ei2 _|_ g2r?2+2J7i+e2i _j_ g2rj2+2r;J+e22 _|_ g2»7i+r?*+»72+ri _|_ g2r?*+»7i+rj2+-ri 
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Figure 3. Non-degenerate type (3,1) incoherently coupled soliton. 



_|_ g2j7i+2»7j-|-»72+r?*+ei2 _|_ g2r?2+2»7*+r?iH-»7|+6»2i _|_ g2r?2+2r?5+»?i+r?^+6»22 _|_ g2(r;i+r;*+??2+r)5)+R4 ('14c) 
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Figure 4. Degenerate type (3,1) coherently coupled soliton. 



and the auxiliary function 

3 3 3 

j=i j=i j=i 

I p''?i+»?2+2''?2+Ai2 _|_ ^ri2+'n'^+2-qi+X2i i „»?2+''?2+2??i+A22 i 2r?i+2r?*+2r?2+Ai 



+ e^ 



2r?i+2r?2+2r?5+A2 , ^2rii+'q*+2ri2+V2+^3 



(Ud) 

In the above equations, rji = ki{t + ikiz), I = 1,2, and the detailed expressions for other 
quantities can be obtained from the Appendix by substituting m = 3. The above (3, 2) 
solution is characterized by eight complex parameters. 

5. Soliton solutions of m— component CCNLS equations 

The bright one- and two-soliton solutions of arbitrary r/i-component system (^ can 
be obtained by generalizing the results of 2-component and 3-component cases. For 
completeness, we give the procedure to write down the multicomponent soliton solutions 
from 771 = 3 case. 



5.1. Bright (m,l) soliton solution 

The (771,1) soliton solution can be obtained from (11) by allowing j to run from 1 to ttt, 
and by redefining kh and e"" as kh = ,^ ^^«. y |q^i T ^^^ e"" 



„(i)\2|2 



respectively. It can also be found that the standard 'sech' type ICS results for the choice 






(i)^2 



and the CCS arises for ^^(«i j 7^ 0. One can also arrive at perfect 
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Figure 5. Non-degenerate type (3,1) coherently coupled soliton. 



(1) (2) 

'sech' type (m,l) CCS for the choice -^^ = -^^ 



(m) 
Am), 
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5.2. Bright (m,2) soliton solution 

In a similar manner, the {m,2) soliton solution takes the form of (14) but with 
j = 1,2,3, ...,m and all the other quantities appearing in the corresponding equations 
are defined in the Appendix. 

6. Soliton collisions in the multicomponent CCNLS equations 

The two-soliton solutions obtained in the preceding sections describe the interaction 
of two solitons in the multicomponent CCNLS system. One can get more insight 
into the collision dynamics and identify interesting collision properties by performing 
the asymptotic analysis of the two-soliton solution. In this section, we analyse the 
2-component and 3-component cases and point out interesting behaviours which are 
not possible in the single component NLS equation. Our analysis can naturally be 
generalized to arbitrary m-component CCNLS system ([3]). To perform the asymptotic 
analysis of the two-soliton solution, we consider the parametric choices kiR, k2R > 
and kij > k2i, without loss of generality. A similar analysis can also be carried out for 
other choices of parameters also. 

6.1. Two-component CCNLS equations 

The 2-CCNLS equations can support both coherently coupled- and incoherently coupled- 
solitons as shown in the preceding sections. Now it is of interest to investigate the 
collision of a CCS with an ICS and also the collision between two CCSs and the collision 
among two ICSs, separately. These three types of soliton collisions in two-component 
CCNLS equations are discussed in this sub-section. 

6.1.1. Collision of coherently coupled soliton with incoherently coupled soliton: First we 
consider the collision of a CCS, say 5*1, with an ICS, say 5*2, in which the former arises for 
the choice («! )^ + ("^i V ¥" and the latter results for the choice (02 )^ + ('^2 V = 0- 
The asymptotic forms of the above solitons 5*1 and 5*2 before interaction {z — )■ —00) 
and after interaction {z — )• +00) can be expressed as below. In the following equations, 
Vi = kit + ikiz), r]in = kmit - 2kiiz), and r]ii = kut + {kfj^ - kfj)z, I = 1,2. 
Before collision (z — t- —00) 
CCS^i {viR^0,V2R^-oo): 

, -, / cos(P-^^)cosh(r77p) + 2 sin(P-^~)sinh(r77p) \ 

where A]- = e ""^'C"\ pj- = ^M_:£, /- = ln(a?^), L^~ = e^^i""^) - 2, and 

"niR = viR+"-f- 

ICS ^2 {ri2R ^ 0, r]iR -^ 00): 

ll = ^sech (r^2R + ^";'" 1 e^^-, J = 1, 2, (156) 
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where A^" = e'^" ~v 2 J. Here and in the following, the superscript and subscript 
appearing in the quantities q, A, P, I, and L represent the soliton number and the 
component, respectively. The sign — (+) appearing in the superscript represents the 
asymptotic form of the soliton before (after) interaction. 
After collision (z — )■ +00) 
CCS Si {viR^0,V2R^oo): 



^'-'^^\ ' Wfe).L.t^ !-"'■ ^■-^■^- <^^=) 



,1+ _ o zii+ I cos(^/^)cos^(^i?) + ^ sin(P/+)sinh(r/+j 

where A]+ = e^- — ), P/+ = ^" , ^S L'+ = e^^'^^V^) - 2, and r]+^ 

ICS ^2 iV2R ^ 0,ViR -^ -00): 




gf = _Z_sech {V2R + ^] e*^^^ j = 1, 2, (15d) 



^i - ^-sech \r]2R + — 

where A^"*" = 02 c~^^. Explicit expressions for various quantities appearing in 
equation (15) can be deduced from the Appendix for m = 2. 

From the above asymptotic expressions (15), we can relate the amplitude (peak 
value of the envelope) of the CCS 5*1 after interaction to that of before interaction as 
A'^-'^ = Tj A^j~, j = 1, 2, where the transition amplitudes Tj's are given by 



1 



^ {kl + k2){h - h)\{a?^22 - aj'^^Ku) + gp^n ' .^^2 (16) 

' V ik,+k*){ki-k*2)^i2\c^?? y ' ^ ' ' 

where f2 = ^ ^ -k ] — ■ ^ careful analysis of equation (I16p shows that the absolute 
values of Ti and T2 will become one only if (al )^ + {a[ )^ = 0, for which one 
can not have CCS. Hence we conclude that in the 2-component CCNLS system ([2]), 
during collision between CCS and an ICS, CCS always experiences an intensity/energy 
switching. In addition to the above amplitude variation, CCS Si experiences an 
amplitude dependent phase shift $1 = ^^^"g^^" = fcT^ ^^ (StT^iItS)- ^^ *^^ 
other hand, the ICS undergoes elastic collision as the amplitude of the ICS ^2 before 
and after collision can be related through the expression A'^'^ = L^_^^w^^~ L A"^" , and 
hence l^^+p = \A^f\'^, j = 1, 2. Also, ICS ^2 undergoes a phase shift $2 = ^^"^2!"^"^" = 

2fei 



^ / "^1- The change in the relative separation distance between these two solitons 
can be written as 

Ati2=tr2-t+2=(l + ^)$l, (17) 

where t]~2 = t2 —t^ and ti2 = t2 —t^ , in which t^ (t^) is the position of soliton Si before 
(after) collision. Here the phase shifts and the relative separation distance between the 
solitons purely depend on fc's but they are independent of a-parameters. 

The main physics behind this collision scenario is that during collision always there 
is an energy/intensity switching in the CCS but the ICS remains unaffected, thereby 



dz 

■ A 

dz 
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leading to energy non-conservation in individual components. However, the total energy 
is conserved. This follows directly from the following expressions 

• / \qi?dt = ^ {qfql-qlqf)dt, (18a) 

J— CO J —oo 

/OO /"OO 

\q2?dt = ^ {qfqf ~ qfql) dt, (186) 

-oo J —oo 

which can be obtained from ([2]) in a straightforward way. In fact, the above relation 
shows that the energy of the individual component, ie. j_ |gjp dt, j = 1,2, is not 
conserved whereas the total energy, that is, J^ (kiP + \q2\'^) dt, is conserved. As a 
consequence of this, the ICS induces significant energy switching in the CCS with an 
amplitude dependent phase shift and reappears elastically after interaction. During 
collision CCS can also switch its profile from single hump to double hump and vice- 
versa. This type of energy switching collision resulting in a dramatic alteration of soliton 
intensity profile is quite different from the shape changing collision of the Manakov 
solitons in HI] and is shown in figure El for illustrative purpose. The parameters are 
chosen as ki = 2.3 + i, ^2 = 2.5 — i, 7 = 2, a\ = 0.75i, a\ = 1.9, Og — 1 + « and 
a2 = I — i- Here, the non-degenerate CCS Si switches its profile from double (single) 
hump to a single (double) hump and exhibits enhancement (suppression) of intensity in 
the gi (^2) component, whereas the degenerate ICS 5*2 undergoes mere elastic collision. 



6.1.2. Collision of coherently coupled solitons: Let us consider the collision of two 
coherently coupled solitons, arising for the choice {a^ )^ -|- (a^- )^ 7^ 0, j = 1,2. Their 
asymptotic forms are found as 
Before collision {z — )■ —00) 
CCS Si. 

^^_ I cos{Pj^)cosh.{r]^j^) + i sm{Pj~)smh.{r]^j^) 
Acosh^ {r]^j^) + L^- 



q]- = 2AY ^^"^^^ ^— V'/i^^ KL ^ I e--'", J = 1,2, (19a) 



*ii'+'. -^11 „i_ <5H'-«- 



where A]- = e^M , pj- = ^iii^, /7 = \n{a?), j = 1,2, L^" = e^R^-"^) -2, and 

CCS ^2: 



cos(Pf )cosh(?72j:j) + i sm{Pj )sinh(?72 
4cosh^(ry2R) + ^^" 



€ = ^^-\ ' .'X-. , rl ^^le--, J = 1,2, (196) 



where A^f = e' "2 , P/" = ^" ^ , L^~ = e^"-(^^) - 2, and r]-j^ 

V2R + ^^T^. 

After collision {z — )■ +00) 

CCS Si . 

1 , 1,1 cos(P/^)cosh(r7i"'p) + i sin(P,^"'")sinh(r7i"*p) , 
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Figure 6. Energy switching collision of a coherently coupled soliton 5*1 with an 
incoherently coupled soliton iS'2 in 2-CCNLS system. 



1+ _ (fcl-fc2)(fci+fc2) 41- 



wnere A- - ^i^._f^,-^^,^^^^,^Aj , t^j 
and 77+fj = r^iij + 



P2/ ^22/ 



P,-, ^^^ 



g022-(^4^) 



L^ 



-R4 — ':22 



CCS ^2: 



,2+ _ o 42+ I cos(P/+)cosh(r72+^) + z sin(P/+)sinh(7y2+^) ^ ^^^^^ , ^ ^ ^ 



2A- 



4cosh2(r/+j) + L2H 



(j) _,+ 



whprP 42+ _ (fci-fc2)(fci+fc2) 42- p2+ _ ^221-^ J I ^ p2- r+ _ ]„(^iJ)\ T 2+ _ JR 



(fci-fc2)(A:i+fc5)""i ' "i 



;i9d) 



4_ 



2 = L , and ?7j^ = rj2R + ^. The other quantities appearing in the above equation 
(19) can be obtained from the Appendix for m = 2. 

From the above equations (llQap -f fTOdj) we find |y4J 



i+|2 



I A 



i-\2 



, hJ 



1,2. 
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Figure 7. Elastic collision of coherently coupled solitons in 2-CCNLS system. 



This clearly indicates that soliton intensities remain same before and after collision. 
Hence the CCSs undergo elastic collision with amplitude dependent phase shift, $i = 

^^'if^ = 41^ ({feSISi) for soliton ^1' a^d $2 = - (If) $1 for soliton 
5*2. Here the change in the relative separation distance between the two CCSs is 
At - ^1 ' '^ 



^12 



1 + f^ ) $1. Such type of elastic collision of the CCSs is depicted in figure [7] 



„(2) 



for ki = 1.5 + i, k2 = 2-i, 7 = 2, al' = 1.7i, a\^' = 1, a'^' = 1i and a)^' = 1.2. 
In figure U\ two CCSs having single hump profiles in qi component and double hump 
profiles in q2 component undergo elastic collision in both the components. 




,j = l,2,(20a) 
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6.1.3. Collision of incoherently coupled solitons: In this subsection, we discuss the 

colhsion dynamics of two ICSs arising for the choices {a, Y + {(^j )^ = 0, j = 1, 2. The 

following expressions are the asymptotic forms of solitons Si and 5*2 before and after 

interaction. 

Before collision {z — )► — oo) 

AYe'"^^' V^ ( sech (r/iR + ^) 

^2-gir,2, j y sech (r/2ij + ^^^) 

where A- = a{ e 2 and A- = e 1 '^ 2 > . 
After collision {z — )• +00) 

'4V^-^Vsech(r^l^ + i^) \ ^._ 12^206) 

Afe^^-' J [ sech {v2r + f ) j '^ " ' '^ ' 

where A- = e 2 ~^ 2 ^^ ^1 T = Og e~^", "T" denotes the transpose of the matrix 
and various other quantities can be obtained from the Appendix for m = 2. In the 
above equation (20) the superscripts 1 and 2 appearing in qj and Aj denote the solitons 
while j denotes the component. 

The amplitudes of the ICSs Si and ^2 after collision can be written from the 
above asymptotic expressions as A]+ = {l^gjgg^" and A^ = Igflf^A^-, 
respectively. One can find the intensities of the solitons before and after collision 
are same, that is lA-^'^P = |v4-~P, ?', / = 1,2. Thus the collision between the ICSs 
arising in the two-component case is always elastic. But Si and 5*2 suffer phase-shifts 

^1 = ^^^^fc^ - 41- ( £;g|g;g ) and $2 = ^^e^ - - (If) $1, respectively, 
as observed from equation (20). From the above phase-shifts one can find the change 
in the relative separation distance as Ati2 = ( 1 + p^ 1 $1. This type of elastic collision 
is exactly similar to the elastic collision of standard NLS solitons and is given in figure 

E for ki = 1.5 + i, k2 = 2-i, 7 = 2, aj^^ = 1 + i, af ^ = I - i, a?^ = 1 + 2i and 

6.2. Three- component CCNLS equations 

Next, we investigate the above three types of soliton collisions in the three-component 
CCNLS equations to examine how they differ from the 2-component case. 

6.2.1. Collision of coherently coupled soliton with incoherently coupled soliton: Here we 
consider the collision between the coherently coupled soliton say Si {{a\ Y + {a\ )^ + 
[af^f ^ 0) and an incoherently coupled soliton ^2 ((«?V + (a?^)^ + ("2^^)^ = 0). The 
asymptotic expressions of the CCS Si and ICS 6*2 before and after interaction can be 
obtained as before and are given below explicitly. 
Before collision {z — t- —00) 
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Figure 8. Elastic collision of incoherently coupled solitons in 2-CCNLS system. 






5 - = 2A 



Acosh^ {r]~j^) + Li- 



e*"", J = 1,2,3, (21a) 



1— „^^ — i pi- "ill 'ji" 



where A- = e 2 , . ^. — 2 ' "j 

ViR = ViR+"-f- 

ICS ^2 (?72i? ^ 0, r^iR ^ cx)): 

q^- = Vsech L« + ^i^^ 1 e^^-, j = 1, 2, 3, 



, /T = \n{a[^^), L^- = e(^i-^) - 2, and 



(21 &) 
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where A^ = e'^" '^ 2 ). 
After collision {z — )• +00) 

1 , 1,1 cos(P/+)cosh(r7it,) + i sin(P/+)sinh(r7i"^„, , 
g;'^ = 24+| ' . X^^^. , l.i ^^i^le^^", J = 1,2,3, (21c) 



4cosh^(r7+j^) + L 



where Aj+ = e^- — ), P/+ = ^" /^S Li+ = e^=^-(^V^) - 2, and 77+^ 



,eu+R2, 



VlR + "-^^■ 



ICS ^2 (^2i?, ^ 0, r/iR -^ -00) 



gf = ^sechh2ii + -^ e^^-, j = 1,2,3, (21d) 



.J = ^-sech ( r/2ij + Y 

where ^4^ = 02 e 2 . 

The amplitudes of the CCS S*! before and after interaction can be related as 
A^-~^ = Tj ^] , j = 1, 2, 3. Here the transition amplitudes T^-'s are defined as 

^ f jkl + h){h - h)\{a?K,, - a^\,,) + a^>n\' \ ' _^2 3 f21e) 

3 
where Vt = ^^^ ^^(aS'^a?^ 

One can check from equation (121 ep that the transition amplitudes become 
unimodular only for the choice {a[ )^ + {a\ )^ + {a\ )^ = 0, for which CCS can not 
exist. Thus the intensities of the colliding solitons before and after collision are always 
different and the CCS Si undergoes energy switching collision with shape alteration 
due to intensity switching among the components and amplitude dependent phase shift 
$1 = ^i^^g^ ^ ^In ( ;:;;gg;| ). But the ICS S, reappears elastically after 
collision with the CCS and its amplitudes in the three components after collision can be 
expressed as A^+ = l^^^gj^j^gj A^f , j = 1,2,3, from which we find \A^+\^ = \A^-\\ 

However 5*2 suffers a phase shift $2 = ^^^2fc"~^" = ~ ( T^j "^i- '^^^ phase-shifts which 
are independent of a-parameters also lead to a change in the relative separation distance 
of solitons Si and 5*2 which is determined as Ati2 = ( 1 + ^^ ] $1. This type of shape 
changing collision scenario is quite different from the shape changing (energy sharing) 
collision of three-component Manakov solitons IHHl] and is shown in figure [9] for the 
parametric choice ki = 1.5 + i, k2 = 2 — i, 7 = 2, ai = 1, a[ = 1.5, al = 2, ^2 
2 + 2, 0^2 = 2 — z, and 0^2 = a/6 i. Here, the CCS 5*1 experiences suppression in its 
intensity after collision and also changes its profile from single hump to double hump 
in gi and qs components but in g2 its intensity gets enhanced and also its single hump 
profile is retained. But the intensity of ICS 5*2 remains intact during collision in all the 
three components. 
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Figure 9. Energy switching collision of a coherently coupled soliton with an 
incoherently coupled soliton in 3-CCNLS system. 
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6.2.2. Collision of coherently coupled solitons: In this subsection, we analyse the 
colhsion behaviour of two coherently coupled solitons arising in the three component 
CCNLS equations for the choice {aff + {aff + {aff ^ 0, j = 1, 2. The asymptotic 
expressions are given below. 
Before collision {z — )■ — oo) 
CCS Si . 

, -, / cos(P/~)cosh(r7rn,) + 2 sin(_P,^~)sinh(?7rR)\ 

(i) — 

where A]- = e '"^'r^" , pj- = ^-^^, Ij = \n{a?), L^- = e^^^-"^) - 2, and 

ViR = ViR+"-f- 
CCS S2 : 

„ „ /cos(P,^~)cosh(r7rp) + i sin(P,^~)sinh(r79p) \ 

g|- = 2A;- ' \ Zl -^ , T2 ^^^ e^^-, J = 1,2,3, (225) 

\ 4cosh [rJ2pi) + L^~ I 

where A^" = e' ""2 , P/" = ^" ^^^" , L^" = e^"-(^^) - 2, and r^s^ = 

_|_ -R4 — £11 

After collision {z — )■ +00) 
CCS Si : 

1, / cos(P,^^)cosh(r7i^,) + z sin(P-^"'")sinh(77i"*R) \ 
,- - 2^- ( ^^^ 4coshV.) + ^-- j ''''^ ' ^ '' '' '' ^^'^^ 

where A^ = ^,^§^A]-, Pj- = ^£# . i^-, L^- = e(^-^^) - 2 . L-, 
and r^+j = r^i^^ + ^^^. 
CCS ^2: 



cos (P^^"*") cosh (r^^) + i sin(P^^"'")sinh(r7^) 
4cosh2(r/+j) + L2H 



gf = 2Af -1-. ;^u.iH,,^;-r.^.iiH^.^,. ;. ,,,^; j ^^^^^ ^ ^, ^ ^^ ^^ ^^ ^^^^^ 



wnere A^- - ^f^^_,^^-^(^,^^_^,^,^A^ ' ^j - 2 - ^j ' S " ^^^"2 J^ ^ - e 2 
2 = L2-, and r/+j = r/a/j + '-f. 

From the expressions for A^-~^ and A^"*" given below equations (122 c\i and (122 (ip . 
respectively, we find that |A'-^P = |^j~P, ^ = 1,2, and j = 1,2,3. This displays the 
elastic nature of collision between two CCSs. Also, the two colliding CCSs Si and S2 
experience phase shifts $1 = ^^^^^ = 4 In ( {^l^gg^g ) and $2 = - (|f) $1, 
respectively. In this case, the change in the relative separation distance between the 
two CCSs is Ati2 = (1 + ^1 $1. Here also the phase shifts and hence the relative 
separation distance of the solitons purely depend on fc's only. This kind of elastic 
collision of the CCSs is shown in figure dUlfor ki = 1.5 + i, k2 = 2 — i, 7 = 2, a\ = 
0.25, af^ = -0.71, af^ = 1.2i, ^^^ = 1, a^ = l.Ai and af = 0.75i. 
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Figure 10. Elastic collision of coherently coupled solitons in 3-CCNLS system. 



6.2.3. Collision of incoherently coupled solitons: The collision dynamics of incoherently 
coupled solitons (say 5*1 and S2) arising for the choice (o;^- )^ + {a\ )^ + (o;^- )^ = 
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0, j = 1,2, discussed here exhibits exciting energy sharing colhsion behaviour which is 
not possible in the two-component CCNLS system. The following expressions are the 
asymptotic forms of the ICSs Si and 5*2 during collision. 
Before collision (z — )• — oo) 

?r V - i (^ ^r^"" ^ " ( «-h ivin + f ) ^,-123 (23a) 

1 .1- (7) _% 1 .2- c(3) ( Ri+R-i -. 

where A^ = a{ e 2 and A^ = e 1 ^ 2 > , 
After collision (z — > 




T 




sech {r]iR + ^^^) 

sech {r]2R + ^) 



,j = 1,2,3,(236) 



where A^^ = e^2 -( 2 ) and A^"'' = cig c ^- ^^ ^he other quantities in equation (23) 
can be obtained from the Appendix for m = 3. 

From the above asymptotic expressions, we arrive at the following expressions 
relating the amplitudes of ICSs before and after collision. 

4+ = T,«4-, A'Y=TJ'^Al, J = 1,2,3, (24a) 

where the transition amplitudes T- and T- of solitons Si and 5*2, respectively, are 
found as 



Ai + 



a. 



i^>n 



,.0) 



K22 



{ki-k2){kl + k2) 






.(2) 



|n|2 1 



^ - ^1 ^2 + ;^77;^ /(A;* - k*){kl + k2 



"2 



Kll 



^r=-^^ ^^ 07 ?c?:::n ^ j = 1,2,3, (24c) 



1_A I "1""" ^ V(^i-^2)(A;i + A;2*) 



2 



3 



where Ai = #H, A2 = #21, and fi = ^.^ ^^(af^a?). 

It can be noticed here that T- and T- , j = 1,2,3, are unimodular only for 

„(!) a<^' q'^^ 

the choice -nj = -fey = -fey and hence the elastic collision occurs only for this 

"2 "2 "2 

choice. Otherwise, the amplitudes of the solitons are different as |T/-'p 7^ 1, / = 1,2, 
j = 1, 2, 3, and this results in an amplitude/intensity redistribution among the two 
solitons split up in the three components. The phase shift experienced by the solitons 
g is $1 = ^^^-f-^i = ^ In ffe^ (1 - Ai A2 + ^^"jl and that of ^2 is found 
to be $2 = — ■^'^i- Here the change in the relative separation distance between the 

two solitons become Ati2 = fl + ^]$i. In contrary to the collision scenario in 
the two component case and also to the two other collision processes discussed in this 
subsection 6.2, here the phase-shift does depend on a-parameters also in addition to /c's. 
Hence the phase shift and the change in the relative separation distance between the 
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colliding solitons can be tuned by altering the polarization parameter a and k suitably. 
This collision scenario is similar to the shape changing collision of solitons involving 
energy sharing among the colliding solitons in the three-component Manakov system 
reported in refs. |H|TT] and will have important applications in the context of optical 
computing and also in achieving multi-state logic [5l[T6]. Figure [TT] shows a typical 
shape changing (energy sharing) collision in which the intensity of solitons 5*1 and 5*2 
is enhanced in gi and g2 components and suppressed in ^3 component after collision. 
Standard elastic collision is depicted in figure [T2j The parameters for figure [TT] are 
h = l.5 + i, k2 = 2-i, 7 = 2, a!^^ = ^2, af^ = V2, af^ = 2i, a!^^ = 



a 



2 = a/6 and ctg = V^- In figurdT2] fci, ^2, 7 are chosen as in figurdTTl and 
a\^^ = 72, aj^ = a/6, aj^^ = VS i. 

This kind of shape changing (energy sharing) collision between the incoherent 
solitons is not at all possible in the two-component CCNLS system as the constraint 
{aj )^ + {a J )^ = 0, j = l,2, restricts them to behave like standard NLS solitons. But 
when we go for three-component the additional freedom involved due to the presence 
of the third component allows the solitons to behave like three-component Manakov 
solitons and results in fascinating energy sharing collisions accompanied by phase-shifts 
which can be tuned by altering both fc's and a-parameters. The shape changing collision 
between coherent soliton and an incoherent soliton in three-component CCNLS is also 
different from the above discussed collision scenario in the sense that in the former case 
energy switching occurs only in a particular soliton (CCS) and the energy of individual 
components is not conserved whereas the total energy is conserved. Also it is not 
possible to achieve elastic collisions during the collision between CCS and ICS. But in 
the present case, energy redistribution occurs in both solitons and also the energy of 
individual components as well as the total energy are conserved. 

6.3. Soliton collisions of ni- component CCNLS equations 

The above analysis of 2-component and 3-component CCNLS system can be extended 
to arbitrary r/i-component CCNLS system. The asymptotic expressions for the three 
different kinds of soliton collisions for the m-component case, with arbitrary m, can 
be obtained from equations (21)-(23) by allowing j and / to run from 1 to m and 
from the expressions given in the Appendix. The collision between the m-component 
coherently coupled solitons again turns out to be standard elastic collision with 
amplitude dependent phase shift and change in relative separation distance. But the 
m-component CCS exhibits intensity switching whenever it undergoes collision with m- 
component ICS and experiences amplitude dependent phase shift which in turn results in 
change in the relative separation distance between the solitons. The transition amplitude 
for the r/i-component CCS is given by 



{kl + k2) {ki-k2)\{a?K22-ai'^ni2 ) + a^^^*nf 
{k, + k*2)ikl-k*2) nl2 \a[' 



T, = { '^'^^'^':r.:'':::^~7 tc ' > ^■=i'2'3,-,-, (25«) 
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Figure 11. Novel shape changing collision of incoherently coupled solitons in 3- 
CCNLS system. 



where Q 



- — ^ ''aj ^2 )• The expressions for phase shift and the change in 



(fci-fca) Z^ 
1=1 



relative separation distance are found to be $i 
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Figure 12. Elastic collision of incoherently coupled solitons in 3-CCNLS system. 
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and Ati2 = (l + ^^ ) $i, respectively. But the incoherently coupled soliton exhibits 
elastic collision in all the m components only with amplitude dependent phase shift. 

An interesting energy sharing collision takes place between the m-component 
incoherently coupled solitons with m > 2 which is not possible in their two- 
component counterparts. In this case both the ICSs undergo shape changing collision 
characterized by intensity redistribution and an amplitude dependent phase shift in all 
the components, which is similar to the soliton collisions of multicomponent Manakov 
solitons [1]. The corresponding transition amplitudes relating the amplitudes of the 
solitons before and after collision can be expressed as 

1 - Ai + ^ 



'^■' .22 



(1) _ V "1 ^2^ ( ih-hm + h) , . _ 1 2 3 m (25b) 



TP-- ''~''''^^^ ( ^'f^:^{ ) > .-M,3...,., (25c) 

.p'*n A \{ki - k2){ki + k^) ' 



A2 + 



Jj) 



where Ai = -4-. — , A2 = -^n — ? and f2 = ,. \. \ 7 («i "2 )• Also, the amplitude 

1=1 

dependent phase shift of the solitons 5*1 and 5*2 can be written as $1 = ^3-R2-Ri ^ 



2/cifl 



In 



2fcifl 
2 /. Of IO|2 






and $2 = — f^*^!, respectively. The above 



k2R. 

transition amplitudes become unimodular and result in elastic collision when the a- 

(1) (2) (m) 

parameters satisfy the specific condition %^ = ^^ = . . . = ^^l_^. 

"2 "2 "2 

7. Conclusion 

In this paper we have considered the integrable ?7i-component CCNLS system ([3]) 
describing simultaneous propagation of m fields in Kerr type nonlinear media and 
obtained the correct bilinear equations by a non-standard bilinearization procedure, 
resulting in more general soliton solutions. The solitons are classified as coherently 
coupled solitons and incoherently coupled solitons depending upon the presence and 
absence of coherent contribution from the co-propagating components/modes. We 
show that one can get coherently coupled or incoherently coupled solitons by tuning 
the polarization parameters suitably. Apart from this, one can also adopt the non- 
standard Hirota's bilinearization method to construct exact soliton solutions of similar 
kind of multicomponent systems which arising in the context of nonlinear optics, spinor 
condensates, etc. Then by considering more general two-soliton solutions we have 
analyzed the collisions among them. Our analysis on their collision dynamics explores 
several interesting collision properties. For the two-component case, we find that the 
CCS(ICS) undergoes elastic collision with CCS (ICS). But during its collision with 
ICS, the CCS experiences energy switching along with amplitude dependent phase shift 
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and change in relative separation distance depending on /^-parameters, leaving the ICS 
unaffected. Our study on the soliton collisions in the three-component CCNLS system 
reveals the fact that the collision dynamics is similar to the two-component case when 
one considers the collision between two CCSs or collision of CCS with ICS. But the 
collision between two ICSs in 3-component CCNLS system displays entirely different 
behaviour from that of two-component CCNLS system. Here both the colliding solitons 
experience energy redistribution and an amplitude dependent phase shift along with 
change in the relative separation distance which depends on both /c's as well as a- 
parameters. This kind of collision scenario is similar to the three-component Manakov 
soliton collisions. Our analysis on the m-component case, with arbitrary m, also shows 
that the collision between CCSs is elastic and always there occurs energy switching 
in the CCS during its collision with ICS. We have also pointed out that the collision 
scenario between ICSs involving energy sharing among the solitons in all components 
for m > 2 and is completely different from m = 2 case, where the collision is mere 
elastic. However the collision between ICSs for m > 2 case, can be made elastic for 
specific choice of a-parameters. This study will find applications in soliton collision 
based optical computing and in optical switches. We believe that our results will also 
have important ramifications in nonlinear optics and in multicomponent Bose-Einstein 
condensates. 
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Appendix 

The various quantities appearing in the sections 3-6 are defined below. 
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Here m, w = 1, 2 and j, / = 1, 2, 3, ..., m. 
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